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A formation scenario for the triple pulsar PSR J0337+1715: breaking a 
binary system inside a common envelope 

Efrat Sabach 1 and Noam Soker 1 

ABSTRACT 

We propose a scenario for the formation of the pulsar with two white dwarfs (WDs) 
triple system PSR J0337+1715. In our scenario a close binary system is tidally and 
frictionally destroyed inside the envelope of a massive star that later goes through 
an accretion induced collapse (AIC) and forms the neutron star (NS). The proposed 
scenario includes a new ingredient of a binary system that breaks-up inside a common 
envelope. We use the binary _c software to calculate the post break-up evolution of 
the system, and show that both low mass stars end as helium WDs. One of the two 
lower mass stars that ends further out, the tertiary star, transfers mass to the ONeMg 
WD remnant of the massive star, and triggers the AIC. The inner low mass main 
sequence star evolves later, induces AIC if the tertiary had not done it already, and 
spins-up the NS to form a millisecond pulsar. This scenario is not extremely sensitive 
to many of the parameters, such as the eccentricity of the tertiary star and the orbital 
separation of the secondary star after the low mass binary system breaks loose inside 
the envelope, and to the initial masses of these stars. The proposed scenario employs 
an efficient envelope removal by jets launched by the compact object immersed in the 
giant envelope, and the newly proposed grazing envelope evolution. 

Subject headings: binaries: close — pulsars: individual (PSR J0337+1715) — stars: 
mass-loss — stars: neutron 


1. INTRODUCTION 


Despite much progress, e.g. review by Ivanova et al. ( 2013(1 . the common envelope (CE) 
phase is one of the least understood processes in stellar binary evolution. The evolution becomes 
much more complicated if a close tertiary body exists in the system, such that the CE evolution 
includes three gravitating bodies. The evolution that formed the Galactic millisecond pulsar (MSP) 
triple system PSR J0337+1715 might have gone through a CE phase involving three stars. 
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PSR J0337+1715 is a triple system that contains a L438 Mq radio MSP with a spin period 
of P = 2.73 ms orbited by two white dwarfs (WDs) (Ransom et al. 2014). The inner WD mass, 
orbital period, and eccentricity are M WD , 2 = 0.197 M 0 , P li2 = 1.63 day, and ei j2 = 6.9 x 10 -4 , 
respectively, while those of the outer WD, hereafter the tertiary star, are M WD ,3 = 0.410M©, 
Pi 2,3 = 327 day, and ei 2i3 = 0.035. The system is coplanar to within 0.01°. Deeper study of the 
inner WD was conducted by Ka plan et ah (2014). 


Taur i s & van den Heuvel (t2014|f propose that the progenitor of PSR J0337+1715 contained a 
massive main sequence (MS) star of mass ~ 10M 0 , orbited by two lower mass MS stars having 
masses and orbital periods of 1.1 M 0 and 835 day, for the inner star, and 1.3M 0 and 4020 day for 
the outer stars. One after the other the two stars entered the bloated envelope of the massive star 
as it became a giant, and ended at very short orbital periods. The neutron star (NS) was formed 
by a core collapse supernova (CCSN) of the primary star. Tauris & van den Heuvel (2014) find 
the range of parameters that can work for their scenario to be about 20%. Later we find that for 
our proposed scenario there is a similar freedom in the parameter space that can work. However, 
we see two problematic steps in this evolution. (1) It is not clear how the two stars ended at short 
orbital periods. If the inner star managed to eject the envelope, as is assumed in the CE evolution, 
then there will be no envelope left to bring the outer tertiary star to a short orbital period. If on 
the other hand there is enough mass left in the primary envelope to swallow the tertiary star and 
bring it closer, then we expect the friction to cause the secondary star to collide with the primary 
core. We avoid this problem by taking the two outer stars to enter the giant’s envelope at the 
same time as a tight binary system. (2) NSs that are bom in CCSNe are bom with a non-zero 
velocity, called a kick velocity. As it is ver y unlikely that the system will stay coplanar after a 
kick, in the scenario of Tau ris & van de n Heu vel (2014) the system gains coplanar geometry at a 
later phase, when the tertiary (outer) star overflows its Roche lobe and transfers some mass to the 
inner binary system, now composed of a NS and a MS star. It is still needed to be shown that such 
a process is capable in bringing a triple system to coplanar geometry. In our proposed scenario 
the NS is formed through accretion induced col lapse (AIC) of an ONeMg WD JNomoto et al.1 
19791: ICanal et al.lll980l: iNomoto & Kondolll991t) . The amount of baryonic mass ejected in this 
process is very small, while the neutrinos are expected to escape in a spherical geometry. No kick 
perpendicular to the orbital plane is expected. The only kick is a small one in the orbital plane due 
to the orbital motion of the WD around the center of mass of the triple system. 

In the standard model MSPs are formed in low mass X-ray binaries (LMXBs) and spin up due 
to ac cretion of mass and angular momentum in a Roche lobe overflow (RLOF) from a co mpanion 
star (IBhattacharva & van den Heuvellll99ll: iTauris & van den Heuvelll2006l: iTaurisI 120151 and ref¬ 
erences therein). In our scenario most likely the tertiary star causes the AIC of the primary star 
that forms a NS, and the secondary star transfers mass to the primary star and causes it to spin up 


and become a MSR Triple star evolution has been proposed for binary pulsars (e.g. Freire et al. 
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2011 : Porteeies Zwart et alJlioi 1 : Piiloo et al.lhoi 2 ). such as the peculiar system PSR J1903+0327 


(Champion et al. 2008). 


Freire et al. 


(1.2011) proposed that two MS stars were spiralling together 


in the envelope of the primary massive progenitor of the pulsar PSR J1903+0327. The inner com¬ 
panion was later destroyed or ejected, leaving behind the binary pulsar. Again, it is not clear that 
one envelope can cause substantial spiralling-in of two well separated stars (the secondary and the 
tertiary). 


The problems we find in the scenario proposed by Tauris & van den Heuvel (2014) motivate 
us to consider a different scenario, also involving a CE phase with three stars, for thejformation of 


PSR J0337+1715. The large fraction of triple systems among binary stars, e.g., Rappaport et al. 


( 2013j) who argue that at least 20% of all close binaries have tertiary companions, further motivates 


us to consider different types of triple systems than previously proposed models for binary and 
triple pulsar systems. In section©we describe our proposed scenario. We employ some well stud- 


evolution when its core mass is ~ 0.2Mm, and more ( 

Refsda 

—- vp 

& Weigert 

1971 

; Webbink et al. 

1983; 

Savoniie 1987; 

Tauris & Savoniie 1999; 

Tauris et al. 2013; 

Ablimit & L 

i 2015; 

Tauris 2015). 


Nevertheless, we also include new ingredients. The new main ingredient, breaking-up a close bi¬ 
nary system inside a CE, is presented in section [3] In section 0 we study the evolution of the 
secondary and tertiary stars to form WDs. Our discussion and summary are in section [5] 


2. THE SCENARIO 


The proposed scenario is^described schematically in Fig©and parameters for one case studied 
here with binary_c ( Izzard et~aill2004 . 2006 . 2009 : section© are listed in Table© The initial triple¬ 
star system is composed from a primary star, the initially more massive MS star, and a tight binary 
system of lower mass MS stars. The outer binary system is tight in the sense that the orbital 
separation of the two outer stars is much smaller than the orbital separation between their center 
of mass and the massive star. The interaction starts when the primary star becomes a giant, and its 
envelope tidally interacts with the tight binary system. 


We divide the evolution of the system into five main stages separated by four phases of strong 
stellar interaction. Here we describe the four interaction phases. The first one is further discussed 
in section© and the last three interaction phases are studied in section© 


(1) Tidal break-up inside a common envelope. The first stage is the most delicate one, as we 
suggest that two processes occur simultaneously: The chaotic process of tidal break-up and the 
poorly understood CE evolution. We start (phase I in Fig© with a triple system composed of three 
MS stars of masses M 1:i w 7.5 — 8.5 M & , M 2A ~ 0.9 — 1.2M 0 , and M 3:i « 1.3 — 1.4, which we 
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tenn the primary, secondary and tertiary, respectively. The two lighter stars form a tight binary 
system, for which we estimate the initial orbital separation for a circular (or low eccentricity) orbit 
to be a 2 , 3 :i ~ 50 — lOOi? 0 . This estimate will have to be checked with a three-body numerical code 
that includes tidal interaction and mass removal by jets from the two stars. This is way beyond the 
scope of this paper. The primary evolves and suffers two expansion episodes, when it has a helium 
core and later when the core is made up mainly from CO. Because of tidal interaction, most likely 
the CE occurs in the first expansion episode (phase II in Fig UK The outer binary system enters the 
envelope of the giant, spirals-in, and removes the hydrogen-rich envelope. 

As the outer binary spirals in, it breaks apart because of tidal interaction. We estimate this 
to occur when the center of mass of the outer binary is at ai j23 ~ a 2)3 ~ 1OO.R 0 from the center 
of the primary. The secondary star, the lighter one, loses angular momentum more efficiently and 
spirals-in to an orbital separation of « 50 — 70i? e from the 2.4M 0 helium core of the primary, 
now having a radius of a few x R 0 . The tertiary star ends further out. This is further discussed in 
section [3] 


(2) AGB interaction. As the primary continues to evolve and forms an ONe core (phase IV 
in Fig{T|), the secondary is as massive as the primary’s envelope and is able to bring it to co¬ 
rotation, enhance the mass loss rate (Tout & Eggleton 1988) and expel the entire envelope of the 
primary (phase V in Fig{T|). Moreover, it is possible that the secondary accretes mass through 
an accretion disc and launches two opposite jets that expel part of the primary He-rich envelope 
(Soker 12014 2015). If a CE prescription is used, like what we will do later, then the removal 
of an envelope by spin-up and jets leads to an effective value of the ac e parameter that is larger 
than unity, acE-eff > 1. At the end of this phase of evolution the primary forms a 1.25 — 1.4M 0 
ONeMg WD, and the triple system has the properties listed in the fourth column of Table Q] for a 
case study. These parameters are justified in the next two sections. In the table, circular orbits are 
assumed for the initial system. 


(3) Tertiary evolution. As the tertiary evolves and becomes an evolved red giant branch (RGB) 
star it interacts with the inner binary system, now composed of an ONeMg WD and a MS star 
(phase VI in FigO}. Two processes take place and lead to the next stage of the triple system that 
we term post tertiary. The first is tidal interaction that leads to the circularization of the orbit, even 
if the orbit is very eccentric. For an eccentricity of ei 2) 3 :p AGB = 0.8 we find the required semimajor 
axis before the evolution of the tertiary (the post primary-AGB stage) to be ai j23:p AGB ~ 55Oi? 0 . 
This implies a periastron distance of « llOi? 0 which leads to a very strong tidal interaction near 
periastron passages and subsequently to a circular orbit. These parameters will bring the final 
result of our simulation even closer to the observed parameters of the system. 


The second process is mass transfer. We find that mass transfer might cause the ONeMg WD 
to overpass the Chandrasekhar mass limit and suffer an AIC (e.g., Tauris et al. 12013c Ablimit & Fi 
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Table 1: Stellar evolution of a case study 


- 

Initial 

post break-up 

post primary-AGB 

post AIC 

final stage 

Observed 

Mi (Mq) 

8 (MS) 

2.4 (He star) 

1.4 (ONeMg WD) 

1.5 (NS) 

1.5 (NS) 

1.4378 

M 2 (M q ) 

1 (MS) 

1 (MS) 

1 (MS) 

1.18 (MS) 

0.2 (He WD) 

0.1975 

M 3 (Mq) 

1.4 (MS) 

1.4 (MS) 

1.4 (MS) 

0.41 (He WD) 

0.41 (He WD) 

0.4101 

dinner [-R©] 

«2,3:i ~ 50 — 100 

Ol,2:BU = 60 

Ul,2:pAGB = 18.3 

ai,2: P Aic = 16.5 

Ol,2:f = 6.3 

6.83 

^outer [Rq\ 

&l,23:i ~ 200 — 400 

«12,3:BU = 214 

fll2,3:pAGB = 270 

ai2,3:pAic = 157 

a 12,3:f = 231 

252.5 

Comments 

e2,3:i = 0 

ei,2:BU = 0.5 

ei,2:pAGB = 0 

ei2,3:pAIC = 0 




ei,23:i = 0 

ei2,3:BU = 0.5 






The masses and orbital separations of the triple system PSR J0337+1715 at the main phases of the evolution, for a 
case study calculated with binary _c (Izzard et al 2004, 2006, 2009). Circular orbits are assumed for the initial setup. 
The table shows an example of one case study with certain initial parameters (column 2; phase I in FigjT} out of the 
allowed parameter range given in the text, e.g., the primary initial mass Mm « 7.5 — 8.5M 0 , secondary initial mass 
M 2 ;i « 0.9 — 1.2Mq, and tertiary initial mass M :i[ \ « 1.3 — 1.4. The post break-up (BU) properties are listed in 
column 3 (phase III in Fig|j}. Note the exchange of inner and outer binary stars from the initial stage to the post 
break-up stage: the initial tight binary is composed of the secondary and tertiary stars, whereas after the break-up of 
the secondary- tertiary binary system inside the common envelope, the primary and secondary stars become the inner 
binary system. The post primary-asymptotic giant branch (pAGB) properties are listed in column 4 (phase V in 
FigjT}, where we find the allowed primary post-AGB mass to range over Mi ;p agb ~ 1.25 — 1.4Mq. During the 
following evolution phase of the tertiary (phase VI in Fig|T} we find with binary jc that the tertiary transfers a mass of 
M a C c:t = O.55M 0 via RLOF to the inner binary. We take 0.44Mq to be accreted by the central binary system taking 
into account mass loss by jets in the accretion disc around the inner binary. Consequently, the primary experiences an 
AIC and forms a NS. We assume its mass is reduced due to energy carried by neutrinos emitted by the cooling NS. 
The post AIC parameters are listed in column 5 (phase VIII in FigjT]). The next stage is the evolution of the secondary 
star, (phase IX in FigjT]) after which the system reaches its final stage (column 6; phase X in FigjT} that corresponds 
with the observed parameters of the system (column 7). 
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2015L for recent papers on the formation of MSPs by AIC). Although there are large uncertainties, 
from binary_c we find that the duration of the RLOF from the tertiary star might allow the AIC 
of Mi. The tertiary star transfers « 0.5M 0 to the central binary system at an average rate of 
pa 1.8 x 10 7 M 0 yr -1 . This will result in an average mass transfer rate o f ^ 1 x 10 7 Mr. : yv~ ] on to 


M. T h is mass transfer is about the rate required to cause an AIC (ITauris et al.ll201 3uAblimit & Li 
20151) . Tauris et al.1 (20131) found that for the case of a giant donor star there is a narrow region 
of donor star masses and initial orbital periods that allows AIC: 0.9 < M donor < 1.1 M@ and 
lOOd < P or b < 300d. Though M 2 is a bit larger than this upper mass limit for M donor the initial 
period and the estimated mass transfer rate found here are in the range required for an AIC. The 
tertiary mass can be fitted to the observed value, as we show in section |U The mass transfer 
increases also the mass of the secondary star and leads to a shrinkage in the orbit of the inner 
binary system (phase VII in FigJT]). In case that the mass transfer from the tertiary star is not 
sufficient for the primary to undergo AIC, only a small amount of mass will accumulate on M\. 
The subsequent evolution of the secondary star, M 2 , will lead the ONeMg WD, M \, to experience 
AIC and become a NS. We do not elaborate on this possibility since we find the AIC via the tertiary 
star to be more likely. 

We take the resulting NS fro m the AIC to have zero kick velocity as is generally agreed upon 
from simulations of an AIC event (IDessart et al.ll2006l:lKitaura et al.ll2006h . Mass loss by neutrinos 
can result in a small kick velocity in the orbital plane of the system due to the orbital motion of 
the WD around the system’s center of mass. The NS symmetrically loses ~ O.15M 0 by neutrinos 
(1/20 of the total mass of the system), with the local NS velocity of ~ 120 km sec -1 (taking the 
system parameters at phase VII in FigJT]). Due to momentum conservation this will result in a small 
recoil of the entire system (« 6 km sec -1 ) in the plane of the system alone, without deviating from 
a coplanar geometry. 

(4) Secondary evolution. The observed mass of the secondary star is 0.1975Mm, and it is a 
He W D. This fits the core of a low mass star during its post-MS evolution (KiDDenhahn & Wciacrt 


1987 


1 99 0), as has been suggested before for low mass binaries (Refsdal & Weigcrt 


1 1971 


Savoniie 


; See section |4] for a discussion on the evolution of the remnant He WD). In our scenario the 


secondary and the NS strongly interact when the secondary evolves as a giant (phase IX in FigJT]). 
According to binary_c the observed orbital separation of ai j2:obs = 6.83 R e dictates the orbital 
separation prior to the AIC phase to be ai i2:p Aic ~ 15.5 — 18.5i? 0 . The mass transfer process will 
occur either through a CE phase (as given by binary_c) or by RLOF from the donor star, M 2 , to 
the central NS, Mi. If only mass transfer through Li is considered, the orbital separation increases 
(e.g., Tauris et al. 20131) . However, the convective envelope and relatively rapid rotation of the 
secondary star result in a magnetically active star that is most likely to posses an extended corona 
with a wind. As the secondary star is already filling its Roche lobe, the corona is now extending 
to the second Lagrangian point, L 2 , and a large fraction of the mass lost from the extended corona 
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is likely to do so from nearby L 2 . As the specific angular momentum of the mass lost through 
L 2 is much larger than that of the binary system, by a factor > 2, the L 2 outflow reduces the 
specific angular m omentum of the inner binary system. This i s will lead to a shr inkage in the 
orbital separation (Livio et al. 1979). From the recent study by Akashi et al. (2015) we estimate 
for the present studied system that by losing & 5% of the total binary mass, ~ 0.1 M 0 , through L 2 
the orbital separation of the inner binary decreases by a factor of ~ 2. This results in a final orbital 
separation of ai, 2:f ~ 6 — 8 R & . 

On the other hand, by evolving the inner binary using binary_c we find that the mass loss from 
M 2 occurs via a CE phase, instead of stable mass transfer via RLOF. The CE phase might be caused 
by the magnetic activity of the secondary star that causes loss of angular momentum from the 
system. To account for the observed (final) orbital separation of the inner binary, ai, 2:0 bs = 6.83_R 0 , 
the CE must be very efficient in removing the envelope. Jets blown by the NS can account for 


«CE-eff > 1- The process might even be a grazing envelope evolution (GEE; Soker 120151; see 
next section) rather than a CE evolution, which implies an even larger value of ctcE-eff (although 
this is not a real CE phase). The tidal interaction and CE evolution explain the observed circular 
orbit. When we use binary jc (section HI), we can justify taking acE-eff > 1- During this phase 
the secondary loses mass, a small fraction of which is accreted by the NS, the orbital separation 
of the inner binary decreases and that of the tertiary increases. Overall we point out that the mass 
transfer can occur in one of the three processes: stable RLOF with some mass loss from the second 
Lagrangian point, an efficient CE phase, or a GEE . The final masses and orbital separations are 
listed in Tabled] 


3. BINARY INSIDE A COMMON ENVELOPE 


The tight binary system (the two lower mass MS stars) tidally breaks-up a pprox imately when 


the o rbital separation of its center of mass from the primary is (IMiller et al.ll2005l: ISesana et al. 
2009h 

3Mi : i(ai, 2 3) 


<*1,23 - 


M, 


2,3:i 


1/3 


°2,3- 


( 1 ) 


The ratio of orbital periods at that epoch is 

-Pi,23 ^ / Q l,23 \ 3,/ ( T7l:i(ai,23) + Af 2;3:i \ ^ ^ ^j 2 / Tfi :i (ai ;23 ) \ // 

p2,3 V a 2,3 / V Tf 2 , 3:i / \-Mi;i(ai, 2 3) + M 2 , 3:i / 


( 2 ) 


where Mi,i(ai, 23 ) is the primary mass inner to the location of the tight binary, ai, 23 , and M 2 , 3:i is 
the combined initial mass of the secondary and tertiary stars. From this condition it is clear that to 
reach tidal break-up the tight binary must spiral-in such that the orbit ai, 23 shrinks faster than the 
orbit of the tight binary system a 2 , 3 . We now discuss the condition for this to occur. 
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The tight binary system loses angular momentum of motion around the primary star inside 
the envelope within a few dynamical times. This implies that the tidal break up occurs on a scale 
not much shorter than the spiralling-in time. Both processes must be considered simultaneously. 
We can safely say that the friction force of the two MS stars in the envelope cannot be neglected. 
If one of the two MS stars is ejected within the envelope it might lose momentum and stay bound, 
even if the star is ejected with an initial velocity above the escape velocity. As we later show, even 
a very high eccentricity of the outer star’s orbit, e ~ 0.8, will lead to the desired outcome. In a 
tidal-break up usually the lower mass companion of the binary is ejected. Inside the CE, however, 
dynamical friction with the envelope must be considered. The outcome might be that the lower 
mass companion of the tight binary is spiralling-in deeper and the heavier companion ends on a 
larger orbit, particularly since the difference in mass between the two MS stars is not large. 


From preliminary simulations (Michaely, E., 2015, private communication) it turns out that 
when the tight binary is immersed in the undisturbed giant envelope the tight binary system be¬ 
comes even harder (its orbit shrinks), and no tidal break-up occurs. The gravitational drag causes 
both MS stars to lose energy, such that the orbital separation between them becomes tighter faster 
than the time on which their distance from the core of the giant star shrinks. In such a case no 
tidal break-up occurs. However, if roughly only the star closer to the giant’s core suffers from 
gravitational drag, tidal break-up can occur. The way to substantially reduce the drag on the outer 
MS of the tight binary system is to remove the envelope outside the orbit of the tight bina ry sys - 
tem. The tight binary system actually experiences a grazing envelope evolution (GEE; Soker 2015) 
rather than a true CE evolution. In the GEE the binary system might be considered to evolve in 
a state of “just entering a CE phase”. The companion star in a regular binary system, or the tight 
binary system in the triple-system case, removes the envelope beyond its orbit by launching jets 
(Soker 2014). The removal of the envelope outside the orbit prevents the formation of a full CE 
phase. This efficient envelope removal increases the effective value of the CE parameter to values 
of «ce > 1 (even though no real CE takes place, but rather a GEE). 


4. FORMING THE PULSAR AND THE TWO WHITE DWARFS 


We use the population nucleosynthesis code binary c of Izzard et al (2004, 2006, 2009) based 
on the Binary Star Evolution (BSE) code of H urle y et al.l (120021) . As binary_c is a stellar population 
nucleosynthesis framework, we use it to calculate the properties of the system at stages beyond the 
tidal break-up of the tight binary, i.e. the mass of the each star, orbital separation, etc. 


To treat the triple system for each phase of strong interaction we simulate a binary system. 
We give a representative case for the scenario, and later show that there is no need here for heavy 
fine-tuning, although the parameter space is narrow. At the post CE break-up phase of the system 
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(phase III in FigJTJ) the primary is a He star of mass M 1B u = 2.4M 0 in an inner binary with the 
secondary star of M 2i bu = 1M 0 . The primary continues to evolve and forms an ONe core in the 
mass range of Mi :p agb ~ 1-25 — 1.4M 0 as it expands along the asymptotic giant branch (AGB; 
phase IV in FigO}. 


Now starts another phase of strong binary interaction. The secondary is massive enough to 
bring the primary envelope to co-rotation, and hence enhances the mass loss rate from the primary, 
up to a factor > 100, e.g.. iTout & Eggletonl (11988b . Another process tha t can aid in removi ng the 
primary He-rich envelope is the launching of jets by the MS companion Jsokerll20f4l ho 1 5l) . This 
stage of interaction with tidally-enhanced mass loss rate and possible jets can not be modelled 
correctly with the binary _c code. We can only mimic these processes by taking the a 0 E parameter 
at this stage to be larger than unity. We here take «ce = 5. From examining other values we 
estimate that to reach the observed parameters of the PSR J0337+1715 system, the inner binary 
separation at the post CE break-up phase should be in the range of ~ 50 — 70 R & . Once the entire 
envelope of the primary is ejected it forms an ONeMg WD, while the secondary star remains a 
MS star (phase V in FigJT]). We term this phase the post primary-AGB phase of the system and we 
estimate the separation of the inner binary at this stage to be ai i2:p AGB ~ 17 — 21 R Q . 


We consider a case study for an example of a possible evolution. We take the immediately 
post primary-AGB system to be composed of an ONeMg WD of mass M 1:pA gb = 1.4M 0 , and 
two MS stars with masses of M 2:pA gb = 1.0M Q and M 3:pA gb = 1.4M 0 . The orbital separation 
of the inner MS star and the ONeMg WD is «i. 2:p agb = 18.3i? 0 , and that of the tertiary star from 
the center of mass of the inner binary is U 3 : P agb = 27Oi? 0 . Circular orbits are assumed for the 
representative case. 


We first follow the evolution of the outer MS star—the tertiary star, as it is more massive 
than the secondary star. For the usage of binary_c we need to treat the inner binary (the ONeMg 
WD + the low mass MS secondary star) as one star when following the evolution of the outer MS 
star. Numerically we take the inner binary system as a compact body of mass M in _ binai . y . pA GB = 
^i, p agb + T/VpAGB = 2.4M 0 , and examine how much mass it accretes from the evolving tertiary 
star (technically we take it as a NS in binary _c). As the tertiary star evolves and fills its Roche lobe, 
it transfers a part of its envelope to the inner binary. To account for the mass loss enhancement 
from the rotation around the massive inner binary ( Tout & Eggletonl 1988 k we take the giant branch 
wind multiplier to be = 10 in binary_c. The tertiary ends this phase of evolution as a He WD of 
mass M 3:f = O.41M 0 ; this will be its final mass. 


With binary jc we find that the tertiary star transfers a mass of M acc:t = O.55M 0 via RLOF 
to the inner binary (phase VI in FigO}- Some fraction, 0.05 — 0.3, of this mass might be blown 
away by jets launched by the accretion discs around the two stars. For that, we take O.44M 0 to be 
accreted by the central binary system. We assume the mass accreted according to the mass of each 
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star; the secondary, M 2 , accretes = O.18M 0 whereas the ONeMg WD, Mi, accretes = 0.26M 0 and 
experiences an AIC to form a NS of mass Mns :p aic ~ 1.5M 0 . As mentioned in section[2l we find 
using binary _c that the duration of the RLOF allows a mass transfer fro m the tertiary to the primary 
at a rate appropriate to cause an AIC (Tauris et al. 2013 : Ablimit & Li 2015). The reduction in the 
mass of Mi is due to energy carried by neutrinos emitted by the cooling NS. Considering the mass 
loss to jets and neutrinos, the orbital separation of the tertiary star from the center of mass of the 
inner binary is somewhat larger than what the binary _c simulation gives. We find the separation at 
the end of the tertiary evolution and AIC stage to be ai 2)3:p Aic = 157 R Q (phase VIII in FigO}- The 
inner binary system is composed now of a NS of mass M 1:pA ic = 1.5M 0 and a MS star of mass 
M 2:p aic = 1-18 M e orbiting each other with an orbital separation of ai j2:p Aic = 16.5/2 0 . 

We next run binary _c for the inner system as the secondary star evolves off the MS and forms 
a He core (phase IX in FigO}. As discussed in section[2l the mass transfer can occur in one of three 
processes: stable RLOF with some mass loss from the second Lagrangian point, an efficient CE 
phase, or a GEE. In the case of a stable RLOF with mass loss from the second Lagrangian point, we 
estimated in section [2] that to reach the final inner binary separation of ai i2: f ~ 6 — 8 R e , a mass of 
~ O.1M 0 should be lost from the second Lagrangian point. On the other hand, evolving the inner 
binary using binary _c results in a CE phase. This CE phase must be very efficient in removing the 
secondary envelope, a process we attribute to jets blown by the NS (Papish et al. 2013 ). Similar to 
the AGB interaction phase, we can only mimic this high efficiency by taking a high value of q 0 e • 
Using binary _c for our case study with «ce = 10 we find at the end of this stage the final masses 
of the NS and the secondary star, now a He WD, to be M NS: f = 1.5M 0 , and M 2:f = 0.2M 0 , 
respectively, and their orbital separation to be ai j2: f = 6.3/7... Due to the mass loss by the inner 
binary system, the orbital separation of the tertiary star increases from ai 2j3:p Aic = 157 R e to 
°i 2 , 3 :f = 231 -Rq. This is the final, and observed, stage of the system (phase X in Fig{l] sixth 
column of tabled]). The masses and orbital separations at the five phases are summarized in Table 

m 


We note that due to some uncertainties in some of the processes, e.g., how much mass is 
lost in the accretion process on to the NS, the final parameters, e.g., final mass of the NS, are not 
certain. Therefore some differences between our derived final parameters and the observed ones 
are not significant. Such is the difference between the observed NS mass and our value in the case 
study, 1.438M 0 and 1.5M 0 , respectively. 

Using the binary_c code for different post primary AGB parameters of the triple system, 
we find that although the parameter space is narrow there is flexibility in setting the post-AGB 
parameters to achieve the observed PSR J0337+1715 system. Namely, there is no need to heavily 
fine tune the post-AGB parameters that can span the following ranges: M 1:p agb ~ 1-25 — 1.4M 0 , 
^ 2 : P agb ~ 0.9 — 1.2M 0 , M 3:p agb ~ 1-3 — 1.4M 0 , ai 2i3:p AGB ~ 260 — 3OO// 0 and ai, 2:p AGB ~ 
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17 — 21i? 0 . We also find that a tertiary post primary-AGB eccentric orbit can be circularized as 
the tertiary becomes a giant and tidally interacts with the inner binary. For example, taking for 
the post primary-AGB tertiary orbit an eccentricity and semimajor axis of ei 2 , 3 : P AGB = 0.8 and 
oi, 23 :pAGB ~ 550/7.,, respectively, leads to the desired post-tertiary parameters. The reason is the 
strong tidal interaction at periastron passages during the RGB phase of the tertiary star. 


The final mass of the secondary star, M 2 = 0.197M©, does not depend much on the post 
primary-AGB separation as well. During the post MS evolution of the secondary star it sub 
stantially expands as a giant when its helium-core mass is ~ 0.2M 0 iRefsdal & Weigert 1971 


Webbink et al.ll 1983c iKippenhahn & Weigertlll990h . and a strong interaction with the NS com¬ 


panion removes the secondary envelope (either in a CE phase or by RLOF, see section [2]). This 
rapid expansion with little increase in core mass has been used before to explain the mass o f 
the He WD com panion of the MSPs ( Savonije 1987 : Tauris & Savonije 1999; Istrate et al. 2014 ). 


Savonije ( 19871) specifically studied a case of binary pulsars that result from mass transfer from 


a (sub)giant to a more massive compact companion. He found that the degenerate helium core 
of the (sub)giant determines the radius and luminosity of the remnant. We point out that in our 
model we concentrate on a scenario where the secondary star had begun evolving after the AIC 
of the NS primary. Though we raise the possibility for AIC from the secondary star, we find the 
AIC earlier through the tertiary star to be more plausible (see section [2]). We use the Modules 
for Experiments in Stellar Astrophysics (MESA), version 5819 (Paxton et al. 2011) to evolve few 
models that are presented in Fig. |2j This figure emphasizes the, well known but crucial for our 
proposed scenario, behaviour described above. Namely, that stars with an initial mass as the sec¬ 
ondary star in our model, fa 0.9 — 1.2M 0 , expand substantially once a helium-core in the mass 
range of fa 0.19 — O.2M 0 is formed. 


5. SUMMARY 


We proposed a scenario f or the triple system PSR J0337+1715 composed of a pulsar orbited 


by tw o white dwarfs (WDs) (IKaplan et al. 


2014: Ransom et al. 2014; Tauris & van den I leuvel 


20141) . Our proposed scenario is described in Fig. HJ and the values for a case study at the different 


evolutionary phases are given in Table Q] 


The triple pulsar system PSR J0337+1715 raises two puzzles. (1) What initial setting can 
lead to a very close inner WD, the secondary, orbiting a pulsar, and a WD, the tertiary, at an 
orbital separation of about the size of a giant star. It is clear that the secondary WD could not have 
evolved through a regular common envelope (CE) phase with the progenitor of the neutron star 
(NS) without the influence of the tertiary star. (2) The triple system is coplanar, indicating that 
no violent event could have taken place, and the system started coplanar. This suggests that the 









































- 12 - 


NS was formed by accretion induced collapse (AIC) rather than by a core collapse supernova. We 
propose that the AIC is caused by mass transfer from the tertiary star during its red giant phase 
(phases VI-VIII in Fig. |T|). Our scenario can also accommodate an AIC induced by mass transfer 
from the secondary star. We do not focus on this scenario since we find the AIC at an earlier stage 
through mass transfer from the tertiary star to be more plausible. 


To solve these puzzles we have introduced to the scenario two new ingredients. The first one 
is the entrance of a relatively close binary system (the tight binary system) into the envelope of the 
primary star during its giant phase. Namely, the two lower-mass stars enter the prima ry env elope 
together, rather than one after the other as in the scenario of|Tauris & van den Heuve l (2014). The 
second new ingredient is the process of grazing envelope evolution (GEE; Soker 2015). 


We assumed that the secondary-tertiary tight binary system evolves in a GEE, such that the 
giant primary envelope outside the center of mass of the tight binary system is removed. The 
spiralling-in process while grazing the giant envelope is a GEE. Namely, the outer star, or a tight 
binary system, is evolving in a ‘just enters a CE phase’. Jets launched by the compact companion, 
or by one or two of the stars in a tight binary, remove the giant envelope outside the orbit. If one is 
to use the «ce parameter, then this efficient envelope removal increases the effective value of the 
CE parameter to values of a CE > 1. In the tight binary system only the star that is closer to the 
center of the primary star suffers the full gas dynamical friction by the envelope. The other star is 
in a region where the envelope gas has been removed already. We consider here an outcome where 
the tight binary spirals-in and tidally breaks-up, such that the lighter secondary star in-spirals and 
forms a binary system with the massive primary core, whereas the heavier tertiary star ends on a 
larger orbit. 


We then used the binary_c numerical code (llzzard et al. 2004120061120091) to further evolve 
the system, first the tertiary star and then the secondary star, as described in section Here as 
well we had to assume that tidal spin-up and/or jets launched by the compact star aid in removing 
the envelope. This implies that we can use a value of «ce > 1- We used values of a CE:IV = 5 
when we evolve the primary star (IV. AGB interaction in Fig. [T]) and a CE:IX = 10 when we 
evolved the secondary star (IX. Secondary evolution in Fig. |T|). It is quite possible that during 
the later CE phases, IV and IX in Fig. [H not only do the jets aid in removing the envelope, but 
the jets are very efficient in doing so and lead to a GEE (Soker 2015[), at least during part of the 
strong interaction phase. On the other hand we do not rule out that the secondary evolution (phase 
IX) occurs through RLOF with mass loss enhancement through the second Lagrange point due to 
the corona of the giant secondary star. Both cases of secondary evolution, either CE or RLOF, 
will result in the shrinkage of the inner binary orbital separation. We also show that no heavy 
fine-tuning is required to achieve the triple system PSR J0337+1715 when one considers tidally 
enhanced mass loss and jet-removal of the envelope. The parameter space to reproduce the system 
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is not large, though. 

Our results have more general implications in showing that a CE with a binary system entering 
the envelope of a giant star can lead to a rich variety of evolved triple systems. The survivability 
of the three stars in such an evolution can be accounted for if the removal of an envelope of a giant 
star with tidal spin-up and jets is considered. 
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Fig. 1.— The evolution of the triple system PSR J0337+1715 from the initial setup, where all three 
stars are on the main sequence (MS), and until the current observations. The initial masses and 
orbital separations represent a case study. There is a range of parameters for which the scenario 
can work. Note that in our notation the initial mass of the tertiary star is larger than that of the 
secondary star: M Xl > M :X l > M 2j . The colors are as follows: grey is for the primary star, 
blue for the secondary star and red for the tertiary star. As each star reaches a more advanced 
evolutionary stage its color darkens. The evolution from phases I to III is assumed here. The 
evolution there after (phases III - X) is calculated with binary jc. Stage IX can occur in one of three 
processes: stable RLOF with some mass loss from the second Lagrangian point, an efficient CE 
phase, or a GEE. We here present the CE possibility since it is the obtained secondary evolution in 
binary jc. 
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Fig. 2.— The evolution of a 0.9 — 1.2M 0 star calculated with MESA (iPaxton et al.l l2 011|) . from 
zero-age main sequence (ZAMS) until the first expansion phase where the star has a helium core, 
as a function of the He core mass. The colors are as follows: blue is for a O.9M 0 star, red is for 
a I.OMq star, green is for a 1.1M 0 star and purple is for a 1.2M 0 star. The thick solid black line 
depicts the radius of the helium core for the case of a l.OM 0 star. The left vertical axis relates to 
the radius of the star (thin solid lines) and the radius of the helium core (thick solid black line), and 
the right vertical axis relates to the He fraction in the core (dashed lines). These plots underline 
the well known behaviour of large stellar expansion in low mass stars of ~ 1M 0 when the helium 
core mass is ~ 0.19 — 0.2M 0 . This has great importance for our newly proposed scenario as this 
behaviour accounts for the mass of the secondary star in the PSR J0337+1715 system. 


core He fraction 























